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Abstract 



We derive a family of ideal (nondissipative) 3D sound-proof fluid models that includes 

j>^' both the Lipps-Hemler anelastic approximation (AA) and the Durran pseudo-incompressible 

'^ i approximation (PIA) . This family of models arises in the Euler-Poincare framework involving 

^ I a constrained Hamilton's principle expressed in the Eulerian fluid description. The derivation 

^^ ' in this framework establishes the following properties of each member of the entire family: 

c/3 . the Kelvin-Noether circulation theorem, conservation of potential vorticity on fluid parcels, 

a Lie-Poisson Hamiltonian formulation possessing conserved Casimirs, a conserved domain 

integrated energy and an associated variational principle satisfied by the equilibrium solutions. 



Having set the stage with the derivations of 3D models using the constrained Hamilton's 
principle, we then derive the corresponding 2D vertical slice models for these sound-proof 
theories. 



> 

in 

^ : 1 Introduction 

o 
en 

The fastest-moving atmosplieric and oceanic waves are the sound waves, whose presence can ad- 
versely affect numerical simulations of atmospheric and ocean circulations, by poisoning the desired 

K> ' low frequency circulations with high frequency oscillations. These high frequency oscillations must 
H ■ somehow be filtered or removed in order to facilitate numerical simulations of the slower motions 

■ ■ ■ that are needed in predictions of weather and climate. An effective way of sound-proofing the nu- 
merical simulations is to replace the exact governing equations with an approximate system that 
does not possess sound waves. Leading candidates for the sound-proofed approximations include 
the anelastic approximation (AA) of [OP62], [LH82] and [Bau96], and the pseudo-incompressible 
approximation (PIA) of [Dur89, Dur08]. Recently, [Kle09] re-examined the range of scales for which 
the AA and PIA sound-proof fluid approximations are valid by using the method of asymptotic 
expansions. 
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Here we derive a family of ideal (nondissipative) 3D sound-proof fluid models from a constrained 
Hamilton's principle in the Euler-Poincare framework. The AA and PIA arise in this family from 
applying different constraints. We then pursue the implications of the framework, particularly 
for conservation of energy and potential vorticity. After analysing the basic properties of this 3D 
family of sound-proof models, we derive their corresponding 2D slice models and analyse their 
properties once again using the Euler-Poincare approach recently developed in [CH13]. These slice 
models assume a |/-independent solution structure but include the y-component of velocity and the 
Coriolis force. Because they assume a constant ^/-gradient of potential temperature, they possess 
solutions that are not solutions of the full three dimensional equations. Slice models have been used 
to study the formation and subsequent evolution of fronts, and are very useful for benchmarking 
numerical schemes since they can be run quickly on a single workstation [Cul07]. 

1.1 Anelastic approximation (A A) 

[Ban96] examined the anelastic approximation for deep fluid convection and proposed an alterna- 
tive form of the anelastic equations. This alternative model combines the results of [DF69] and 
[LH82] to produce a hybrid theory that (1) conserves the domain integrated energy; (2) preserves 
potential vorticity on fluid parcels; and (3) accurately represents the acoustic adjustment process in 
Lamb's problem, cf. [Ban95]. These equations for a dry anelastic compressible fluid (atmosphere) 
rotating at angular frequency f2 under constant vertical gravitational acceleration gz comprise the 
following constrained dynamical system, 

^ + 20xu = -v(i^)+fs, (1.1) 

at \poJ Oq 

-"V^') = 0. (1.2) 

dt ' ^ ' 

with the density-weighted incompressibility constraint 

V-(poWu) = 0, (1.3) 

and two diagnostic relations for the thermodynamic variables p' and T', 

^' P' P' ^ P' P' T' 

— = —^ — , and — = — H . (1.4) 

9q pogHp po po po To 

In these equations, the fluid velocity is denoted u, the advective time derivative is written as 
d/dt = d/dt + u ■ V, and the state variables for this anelastic motion are: pressure p, density p, 
potential temperature 9, and temperature T. Each of these state variables consists of the sum of 
the base state (with subscript 0) and a dynamic contribution denoted with a prime, as in 

e{x,y,z,t) = eo{z)+e'{x,y,z,t), 
D{x,y,z,t) = po{z)+ p'{x,y,z,t), etc., 

where 6' (resp. p') is the dynamic contribution to the potential temperature (resp. mass density) 
field. The base state satisfies 

-T- = -9pQ-, Po = PoRTo, C^9o— = -g, (1.5) 
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where ttq = Tq/6q. The constants R and C^ are the ideal gas constant R and the specific heat at 
constant pressure C„ for dry air. The scale height of the base density state is 

l/H, = ~poiz)-'dpo/dz. (1.6) 

Given the base state functions satisfying relations (1.5), as well as the velocity u and the dynamic 
contributions p' and 6' at any time, the system is completed by the constraint in (1.3), whose 
preservation determines p', and the thermodynamic diagnostic relations in (1.4). The distinctions 
between the anelastic approximation (AA) and the traditional models [DF69], [LH82] are discussed 
in detail by [Ban96] and [Kle09]. The important point for us here is that the model equations (1.1) 
- (1.6) agree in both formulations, that of [LH82] and that of [Ban96]. 

1.2 Pseudo- incompressible approximation (PI A) 

The dynamical equations of the pseudo-incompressible approximation (PIA) are given by [Dur89] 
as a motion equation and an advection law. These are, respectively, 

^ + 2f2 X u = - (^0 + 0')Vp' + ^z , (1.7) 

a„d ^J^2^ = 0. (1.8) 

dt 

Note that the pressure force in the PIA motion equation (1.7) differs from that in the AA motion 
equation (1.1). The weighted incompressibility constraint for PIA is 

V-{poiz)eoiz)n) =0, (1.9) 

which also differs from the corresponding density-weighted incompressibility constraint (1.3) for 
AA. However, both the entropy (i.e., heat and buoyancy) advection law (1.8) and the two diagnostic 
relations (1.4) for the thermodynamic variables are the same for the two models AA and PIA. 

The AA and PIA equations differ in the pressure forces in their motion equations, (1.1) versus 
(1.7), and in their weighted incompressibility constraints, (1.3) versus (1.9). As we shall see, these 
differences are related. Namely, the difference in pressure forces arises from a difference in the 
constraints enforced by the pressure as a Lagrange multiplier, in the Hamilton's principles for the 
two theories. 



1.3 Outlook 

In this paper, we show that both sets of motion equations (1.1) for AA and (1.7) for PIA are 
members of a single family of Euler-Poincare equations obtained from a constrained Hamilton's 
principle expressed in the Eulerian fluid description. Two sets of constraints arise from the under- 
lying Lagrangian fluid flow. One set of constraints affects the variations of velocity, density and 
potential temperature. The other set of constraints imposes sound-prooflng by using the pressure 
as a Lagrange multiplier to enforce conditions on the density and potential temperature. The Euler- 
Poincare formulation of this family of sound-proofed fluid equations provides their Kelvin- Noether 
circulation theorems which lead to the conservation of potential vorticity on fluid parcels. Domain- 
integrated energy is also conserved for the entire sound-proofed family of equations. This result can 
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be demonstrated by applying a Legendre transformation to obtain a Lie-Poisson bracket formula- 
tion of the sound-proofed family of equations. In turn, the Lie-Poisson Hamiltonian formulation 
reveals a set of conserved integral quantities known as Casimirs. In the underlying Lagrangian 
fluid flow, the Casimirs generate motion of the fluid along streamlines of the steady solutions of the 
sound-proofed family of equations. Thus, the Casimirs provide a constrained-energy variational 
principle for obtaining equilibrium solutions of sound-proofed fluid dynamics. They also form the 
basis for determining the Lyapunov stability conditions for these equilibria, following the approach 
applied to the Euler-Boussinesq equations in [AHMR86]. (The Euler-Boussinesq equations form 
a special case of the sound-proofed models in which the base state is constant.) In general, these 
conservation laws are very important for the analysis of geophysical flows and the design of numer- 
ical schemes for weather prediction and climate modelling. For example, in [Ciil07], conservation 
laws were emphasised in analysing the ability of numerical schemes to reproduce balanced model 
solutions in asymptotic limits. See [ThuOS] for a review of the relevance of conservation laws to the 
design of dynamical cores of weather and climate models. The Euler-Poincare formulation allows 
us to produce hierarchies of models that all share the same conservation properties, from fully 
compressible deep atmosphere models to semi-geostrophic balanced models, for example. Here we 
show that the Euler-Poincare formulation also unifles the family of sound-proofed equations. 

Previous related work. In previous related work on Hamiltonian structures of soundproof 
models, a two-dimensional study of the Hamiltonian structure of the Lipps-Hemler AA model was 
presented by [SS92], who also studied wave-activity conservation laws in the two-dimensional case. 
A canonical Hamiltonian formulation of the 3D Lipps-Hemler AA model in the Lagrangian fluid 
description was given in Appendix A of [Bcr95]. As one would expect, the Lie-Poisson Hamiltonian 
formulation for the Eulerian fluid description of these equations is equivalent to the canonical 
formulation of [Ber95] in the Lagrangian fluid description. This equivalence is guaranteed by the 
Euler-Poincare theorem [HMR98]. This theorem was developed for applications in geophysical 
fluid dynamics (GFD) in [HMR02] and specialised to slice geometries in [CH13]. It states that the 
following four dynamical perspectives of fluid mechanics are equivalent: (1) Hamilton's principle 
for the Lagrangian fluid description; (2) the Euler-Lagrange equations in the Lagrangian fluid 
description; (3) Hamilton's principle for the Eulerian fluid description with certain constrained 
variations; and (4) the Euler-Poincare equations in the Eulerian fluid description. 

The methods of this paper are based on the Euler-Poincare theory of reduction by symmetry 
of variational principles introduced in [HMR98]. However, in this paper we shall just quote the 
main results of this theory and apply them to derive a family of sound-proofed models. 

Organization of the Paper. In §2 we recall from [HMR98, HMR02] the results of the Euler- 
Poincare theorem for Lagrangians in continuum mechanics depending on advected quantities (such 
as density and potential temperature), along with their associated Kelvin-Noether theorem and 
Lie-Poisson Hamiltonian formulations. These results establish the mathematical framework into 
which we place the dynamical equations for the AA and PIA models in §3. In §4 we introduce 
soundproof slice models that preserve the fluid properties guaranteed by the Euler-Poincare frame- 
work. Summary conclusions of the main results of the paper and some of the outstanding problems 
it raises are discussed in 55. 
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2 The Euler— Poincare Theorem in GFD 

We begin by recalling from [HMR98, HMR02] the statements of the Euler-Poincare equations and 
their associated Kelvin-Noether theorem in the context of continuum mechanics and approximate 
models in GFD. 

The Euler-Poincare equations for a GFD reduced Lagrangian £[u,D,6] (usually the domain 
integrated kinetic energy minus the domain integrated potential energy) involve the Eulerian fluid 
velocity vector u, the scalar potential temperature 6 and the mass density D as functions of three 
dimensional space with coordinates x and time t. In vector notation, these equations are expressed 
as [HMR98, HMR02], 

or, equivalently, in curl form as, 

= -— I ) — u X curlf 

dt\D6uJ \DSu 

where Sl/SD is the variational derivative of / with respect to D, defined by 



/ 



Si ^^^ ,. i[u,D + e6D,e]-i[u,D, 

- — oU ax = lim — 

5D e^o e 



for all density perturbations 6D. Similar definitions hold for Sl/66 and 6l/du (the latter of which 
is the linear momentum, a vector quantity). The Euler-Poincare system is completed by including 
the auxiliary equations for advection of the total potential temperature 6, 

an 

— + u-V9 = 0, (2.3) 

and the continuity equation for the mass density D, 

^ + V-(Du) = 0. (2.4) 

Specific models are then obtained (for example, the sound-proof models discussed in this paper) 
by choosing a particular form of the Lagrangian i. 

Kelvin-Noether circulation theorem. The Euler-Poincare motion equation in either form 
(2.1) or (2.2) results in the Kelvin-Noether circulation theorem, 

d I ^ Si ^ /I 5i^^ , 

— (* — — •rfx=-(* — — Ve-dx, 2.5 

where the curve 7t(u) moves with the fluid velocity u. Then, by Stokes' theorem, the Euler- 
Poincare equations generate circulation of the quantity D~^5i/5\\ whenever the gradients V^ and 
V{D~^6i/S6) are not coUinear. 
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Potential vorticity conservation laws. Taking the curl of equation (2.2) and using advection 
of the potential temperature 6 and the continuity equation for the density D yields conservation 
of potential vorticity PV on fluid parcels, as expressed by 

^ + u-Vg = 0, (2.6) 

where q = — V6' ■ curl I — -^ 1 . (2.7) 

D \DouJ 

Consequently, the following domain integrated quantities are conserved, for any smooth function 



C$= f d^x D<^{9,q), V$. (2.8) 



Legendre transform, energy conservation and Hamiltonian formulation. The absence 
of explicit time dependence in the Lagrangian i [u, D, 6] gives the conserved domain integrated 
energy, via Noether's theorem for time translation invariance. This energy is easily calculated 
using the Legendre transform to be 

/r /J 
d^x (u-—)-i[u,D,9]. (2.9) 

When the Legendre transform is completed to express E [u, D,6] as if [m, D, 6] with 

m = 5i/6u and 6H/6m = u, (2.10) 

the Euler-Poincare system (2.1)-(2.4) may be expressed in Hamiltonian form as 

^ = {fi,H}, with fiE[ni,D,9], (2.11) 

and the Lie-Poisson bracket is given in Euclidean component form by 

{F,H}[in,D,9] = 



J '^^'' [j2-§^. 5^ (o'i"^* + "^j-^^ ) 



SH 

Sm,j 






(2.12) 



J ■' j ■' ^ 



The conserved quantities C$ in (2.8) are then understood in the Lie-Poisson Hamiltonian 
formulation (2.11) - (2.12) of the Euler-Poincare system (2.1) - (2.4) as Casimirs that commute 
under the Lie-Poisson bracket (2.12) with any functional of (m, D,6'). The Casimirs also result 
via Noether's theorem from symmetry of the Hamilton's principle for the Euler-Poincare system 
under the particle relabeling transformations that leave invariant the Lagrangian i[u,D,6]. For 
full mathematical details, consult [MR95, HMR98, HMR02, CH12]. 
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The Euler— Poincare framework. The four properties (2.5)-(2.9) and the Lie-Poisson Hamil- 
tonian formulation (2.11) - (2.12) of the Euler-Poincare equation (2.1) and its auxihary equations 
(2.3) and (2.4) are all desirable elements of approximate models for applications in geophysical 
fluid dynamics expressed in the variables {u,D,6). Thus, the Euler-Poincare theory offers a uni- 
fied framework in which to derive approximate GFD models that possess these properties: the 
Kelvin-Noether circulation theorem, conservation of potential vorticity on fluid parcels, and the 
Lie-Poisson Hamiltonian formulation with its associated conserved Casimirs and conserved domain 
integrated energy. Previous work [HMR98, HMR02] has shown that many useful GFD approxima- 
tions may be formulated as Euler-Poincare equations, whose shared properties thus follow from 
this underlying common framework. The aim of the next section of this paper is to cast the sound- 
proof motion equations, including those for AA (1.1) and for PIA (1.7), into the Euler-Poincare 
framework. 



3 Hamilton's principle for sound-proof motion equations 

3.1 The sound-proof Lagrangians 

In the Eulerian fluid representation, we consider Hamilton's principle SS = for fluid motion in 
a three dimensional domain with action functional S = j dti and Lagrangian i[u, D,6]. For the 
sound-proof models, we propose a Lagrangian given by the total kinetic energy in the rotating 
frame, minus the total potential energy including the thermodynamic energy, plus the sound-proof 
constraint imposed by the pressure p' as a Lagrange multiplier. 



t;sp 



D[ -\u\^ + u ■ R{x) - gz - C^iTo{z) e 



+ p'{po{z)o{eo{z))-De{e) 



d^x. (3.1) 



SP constraint 



Here D denotes the mass density, u is the Eulerian fluid velocity, R is a vector fleld whose curl is 
2f2 (twice the local angular rotation vector), 6 is the total potential temperature, ^o(^)) 7'"o(^) and 
po{z) are reference proflles and O is an arbitrary smooth function still to be chosen. 

The constraint imposed by the pressure. In the Lagrangian (3.1) the pressure p' plays the 
role of a Lagrange multiplier in Hamilton's principle that enforces a constraint that relates the 
mass density D and potential temperature 6. To understand the meaning of this constraint, one 
may combine (2.3) and (2.4) into 

?a,^.^O0ie).)^O, (3,2) 

where 6* is a smooth function. Imposing the sound-prooflng constraint D0{6) = pq{z)0{6q{z)) 
then yields 

V-{po{z)0{eo{z))u)=0. (3.3) 

For 0{6) = 1 — a + a6, one recovers PIA flows for a = 1 and AA flows for a = 0. For a = and 
Pq{z) = const, one recovers divergence-free flows. 
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Time-dependent reference states. Following [Dur08] one may allow reference states with 
prescribed time-dependence by replacing po{z) by p{x,y, z,t), 9o{z) by 9{x,y,z,t) and tiq^z) by 
7T{x,y, z,t). Imposing the sound-proofing constraint D0{6) = pO{6) then yields 

v-{pe{e)u) = -dt{pO{e)), (3.4) 

in which the right-hand side is a prescribed function of (x, y, z, t). The pressure is then determined 
from the sound-proof motion equation via an elliptic equation involving 

V ■ {pO{9)dtu^ = -dl [pO{e)) - V ■ (udt{i50{e)) 

The resulting sound-proof fiuid equations keep their form, as long as the time- dependent back- 
ground state remains in hydrostatic balance. However, the time-dependence in p, 9 and tt prevents 
conservation of energy, although potential vorticity is still conserved on fiuid parcels. 

Variational derivatives of the Lagrangian fsp- The Lagrangian £sp in (3.1) possesses the 
following variational derivatives at fixed x and t. 

u + R X 



D 5vi 

up + u ■ R(x) — gz 



SisP 1, ,2 



5D 2' 

-C,7io{z)9-p'Oi9), (3.5) 

1 54p 



D 69 

Sisp 



Sp' 



-C,noiz)-p'0'i9), 
Po{z)Oi9oiz))-DOi9). 



One obtains the corresponding motion equations upon substitution of these variational deriva- 
tives into the Euler-Poincare formula (2.1). 

3.2 The motion equation for SP 

From the Euclidean component formula (2.2) for Hamilton's principles of this type, we find the 
fluid motion equation for the SP model in three dimensions, 

^ - u X curlR + 0{9)Vp' +(g + CJ^] z = , (3.6) 

at \ dz J 

where curlR = 2i7(x) is the Coriolis parameter (i.e., twice the local angular rotation frequency). 
We use equation (1.5) to rewrite the last term in parentheses as 

, + C„e- = g[l--y-g^^. (3.7) 

Consequently, we obtain the SP motion equation, namely, 

^-ux2f2(x) + 6'(^)Vp'-^z = 0, (3.8) 

dt 9q 

as the Euler-Poincare equation for the SP Lagrangian (3.1). 
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Solving for the SP pressure. The SP contraint with arbitrary smooth function 

Poiz)Oieo{z)) = DOi9) (3.9) 

is obtained from stationarity of the Lagrangian (3.1) with respect to variations in p'. Upon sub- 
stituting this constraint into the continuity equation (3.2), we find the SP divergence condition in 
equation (3.3). Preservation of this condition determines the pressure p' analytically, by solving 
the elliptic equation obtained by taking the divergence of the SP motion equation (3.8) after first 
multiplying it by the product Pq{z)0{9o{z)). 

The boundary condition for the resulting elliptic equation is obtained from the normal compo- 
nent of the SP motion equation (3.8) evaluated on the boundary and using the boundary condition 
for the velocity, e.g., that it has no normal component at the boundary, which yields a Neumann 
boundary condition for obtaining the pressure. 

3.3 The Kelvin— Noether theorem for SP equations 

From equation (2.5), the Kelvin-Noether circulation theorem corresponding to the SP fiuid motion 
equation (3.8) in three dimensions is, 

— (f> (u + R) ■ dx 



/ {0{9) Vp' - O^Wz) ■ t/x , 

[he X {o'{e)Vp - |-Vz) ■ s 

JJs ^0 



(3.10) 



where the curve 7i(u) moves with the sound-proof fiuid velocity u. By Stokes' theorem, these 
sound-proof equations generate circulation of (u + R) around 7i(u) whenever the vectors V6', 
Vp' and V-2 are not collinear. Using advection of 9 and the SP continuity equation, one finds 
conservation of potential vorticity gsp on fiuid parcels, cf. equation (2.6), 

^ + u-Vgsp = 0, (3.11) 

where gsp = 7^^^ ' curl (u + R) . 

Consequently, the following domain integrated quantities are conserved, for any smooth integrable 
function $, cf. equation (2.8), 

C^= [ D<^{9,qsp)d^x, for all $. (3.12) 
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3.4 Energy conservation, Lie-Poisson Hamiltonian formulation and 
nonlinear Lyapunov stability analysis 

One uses the Legendre transform (2.9) of the Lagrangian (3.1) to find the corresponding Hamilto- 



nian, 



- p' (^poiz)ei9o{z)) - D 0(^)) d'x . (3.13) 

SP constraint 

The Lie-Poisson bracket (2.12) now generates the SP motion equation (3.8), as well as the 
auxiliary equations (2.3)-(2.4) and the SP constraint (3.9) from the Hamiltonian Hsp according 
to equation (2.11). The conserved energy is obtained by evaluating the Hamiltonian Hsp on the 
SP constraint, as 

Esp= f (^\u\^ + gz + C,7ro{z)9\Dd''x. (3.14) 

This is the sum of the kinetic energy, gravitational potential energy and thermodynamic internal 
energy. That is, the SP constraint (3.9) does not interfere with conservation of the usual total 
energy. 

Variational principle for SP equilibrium solutions. The equilibrium solutions of the SP 
equations occur at critical points of the sum if$, where 

H^ = Hsp + C$ , (3.15) 

and C$ represents a family of Casimirs given by 

C^= [ D<!>{9,q)d^x, V<l>. (3.16) 



Here $ is an arbitrary smooth function and q is defined upon using the definition of the momentum 
density m in (2.10) as, 

q = —V6-cml{m/D). 

Thus, the Casimir conservation laws for our Lie-Poisson Hamiltonian formulation of the three- 
dimensional anelastic equations in the Eulerian fiuid description provide a constrained-energy 
variational principle for the equilibrium solutions of the dynamical SP equations and form a basis 
for determining their Lyapunov stability conditions, as done for the Euler-Boussinesq equations 
in [AHMR86]. The AA and PIA equations comprise special cases of the SP equations when 
0{9) = 1 — a + aO. In turn, the Euler-Boussinesq equations form a special case of the AA and 
PIA equations in which the base state is constant. Consequently, the analyses of the nonlinear 
Lyapunov stability conditions for the equilibrium solutions of the three-dimensional SP equations 
and their special cases the AA and PIA equations all follow a similar procedure to that performed 
in [AHMR86], which will produce a similar result, modified according to the non-constant base 
state. 
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4 Sound-proof vertical slice models 

Slice models are used to describe the formation of fronts in the atmosphere and ocean, modelling 
the situation in which there is a strong horizontal temperature gradient which maintains a vertical 
shear flow in the direction tangential to the gradient, through geostrophic balance. On the f -plane, 
this base flow can be modelled with a three- dimensional flow with constant temperature gradient 
in the y-direction, and velocity pointing in the x-direction with a linear shear in the z-direction. 
This base flow is unstable to y-independent perturbations to all three components of velocity and 
temperature, rapidly leading to y-independent front formation. The y-component of velocity is 
coupled to the other variables through the constant y-gradient of temperature. In the case of the 
incompressible Euler- Boussinesq equations, solutions of these equations are also solutions of the 
full three-dimensional equations. In other cases, solutions of the three-dimensional equations are 
not recovered but the slice models provide very useful tools for analysing numerical methods for 
the atmosphere since they can be compared with semigeostrophic solutions [Cul07]. 

[CH13] considered variational formulations for vertical slice models in which the prognostic 
variables are independent of y, with the exception of potential temperature, which has a time- 
independent constant y- derivative, i.e., 

e{x, y, z, t) = es{x, z, t) + {y- yo)s, (4.1) 

where s is a constant. Here, we adopt the notation u = {u,w), and treat v separately, where u, 
V and w are the x- y- and 2;-components of velocity, respectively. We also write V = {dx,dz). 
Consequently, the ^/-independent mass density D{x,z,t) satisfies 

dtD + V- (uD) = , (4.2) 

and three-dimensional scalar tracer equation (2.3) becomes a dynamic equation for 6s{x, z, t) which 
satisfies, 

dtOs + u-Ves + vs = 0. (4.3) 

For a shce Lagrangian l[u,v,D,9s], [CH13] showed that these conditions lead to the slice 
Euler- Poincare equations (written here in curl form), 

d I 5(. ( \ bl 

7— -;- -: U X curl -- -^ 

dtDbw. \ Bbw. 

where in the shce notation we write txxcurl v = {—w{dzVi—dxVs),u{dzVi—dxV3)). The system (4.4- 
4.5) is completed by including the advection equations (4.2) and (4.3) for D and 9s, respectively. 
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Kelvin-Noether theorem. [CH13] showed that the equations (4.4-4.5) can be combined to 
obtain 

--m + uxcurl(^-mJ+V^n.m--J = 0, 

where 

— SI 51 „, 

m = s- ^^Os- 

ou ov 

This leads to the Kelvin-Noether conservation law for circulation obeyed by slice models in the 

form 

where c(u) is any closed material loop that moves with the slice fluid velocity u in the vertical 
slice plane {i.e., not in the ^/-direction). 

Potential vorticity Equation (4.6) implies that potential vorticity q is conserved along flow 
lines of the fluid velocity u, 

dtq + u-Vq = 0, (4.7) 

with potential vorticity 

g = V^m = sV^ . (1£) - V^£ ■ V^^, (4.8) 

which is in fact the usual Ertel PV written in the slice variable notation. 

Energy conservation As for the three-dimensional case, the absence of explicit time dependence 
in the slice Lagrangian i[u,v,D,6s] gives the conserved domain integrated energy, via 
Noether's theorem for time translation invariance. This energy is again easily calculated using the 
Legendre transform to be 



(u- — + v-—)dxdz 

-£[u,v,D,es]. 



When the Legendre transform is completed to express E [u, v, D, ds] as if [m, n, D, 9s ] with m = 
S£/6u, n = bljbv, 5H/5\w = u, and 6H/6n = v, the Euler-Poincare system (4.2-4.5) may be 
expressed in Hamiltonian form as 

^ = {fi,H}, with fie[m,D,9s], (4.9) 
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and the Lie-Poisson bracket is given by 

{F,H}[ni,D,9] = 

((— v^— -1^— V^— ^d 
\\5rn, / 5m \5m J 5m J 

, { 5F ^5H 5H ^5F\ , , 

- I n{^- V- —- ■ V^ dxd^ 

\07Ti 5n 5m 5n J 

5F ^^ 5F \ 5H ^ ^ 

5m 5n J 56s 

5H ^^ 5H \ 5F ^ ^ 

5m 5n J 56s 

4.1 Hamilton's principle for sound-proof slice models 

In the present notation, the Lagrangian for the Shced Sound-proof Approximation (SSA) in Eule- 
rian (a;, z) coordinates is, 

Issa[^.v,D,6s] = / -{\u\^ + v^)+fDvx-gDz 

-DC,noiz)6s (4.10) 

+ p'(poiz)Oi6'si^))-DOi6s)) dxdz, 



SSA constraint 

As before, the constraints are imposed by the pressure as a Lagrange multipher. 
The variational derivatives of the SSA Lagrangian are, 

u, -f^—, — = v + fx, 



D 5u ' D 5v 

^^SSA 1 /I |2 I 2\ . r 

— — — = - u + V ) + fvx 
5D 2 ^' ' ' ■' (4.11) 

-gz-C,no{z)6s~p'0{6s), 

1 ^^SSA 



c,Mz)-p'e'{6s). 



D 56s 
For the shced sound-proof approximation SSA, the Euler-Poincare approach yields the following 
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equations on the slice semidirect product with advected density D and scalar 9s'- 

T—U + u ■ Vu — fvx 

ot 



+ (9 + c.^Os) z + oies)W = 0, 



dv ^ „ „ 

- + u.Vv + x.uf ^4_^2) 

-{C,no{z)+p'0'{es))s = 0, 
f + V.(.D)^0 

ot 

For 0{6) = 1 — a + a6, we have SPIA flows for a = 1 and SAA flows for a = 0. For a = 
and po{z) = const, one recovers sliced incompressible flows. Note that the solutions of the slice 
equations for SSA models are not solutions of the full three-dimensional equations. 

One obtains the expected Kelvin circulation conservation law (4.6), written for the SSA models 
as 

i. / {su-{v + fx)Ves)-dx = 0, 

"'^ Jc{u) 

which leads to potential vorticity conservation 

dtq + u- Vg = 

with potential vorticity 

g = sV^ ■ u - V^(w + fx) ■ VOs, 

This is the same circulation theorem as for the Euler-Boussinesq slice model in [CHI 3]. 

The Euler-Poincare slice equations are Hamiltonian, with conserved energy 

H= I ^{\m\^ + {n-Dfxf) 

+ gDz + DC^'KQ{z)es d^x , 
and Casimir conservation laws, 



C$= / D$(g) dV, (4.13) 

for an arbitrary smooth function $. Furthermore, critical points of the sum 

Hc = H + C^ (4.14) 

are equilibrium solutions. This is a general feature of Lie-Poisson Hamiltonian theories, see, e.g., 
[HMRW85], and in this case is the basis for studying Lyapunov stability for critical-point equilibria 
of the slice models. However, this feature will not be pursued further in the present paper. 
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5 Conclusions 

In this paper we derived in the Euler-Poincare variational framework a family of 3D soundproof 
models that contains the Lipps-Hemler anelastic approximation (AA) and the Durran pseudo- 
incompressible approximation (PIA) as special cases. In this family of 3D soundproof models, 
the pressure arises as a Lagrange multiplier enforcing a constraint on the density, just as in the 
standard incompressible case. The models in the sound-proof family only differ in how the con- 
straint is enforced. We have explained how these constraint forces relate to the standard forms 
of the equations. Having found Lagrangians that lead to these models in Euler-Poincare form, 
we immediately deduce that the models have a conserved energy, inherit a Lie-Poisson structure, 
and have a Kelvin-Noether circulation theorem and corresponding conserved potential vorticity. 
This formulation opens the way for further approximations in the Lagrangian, such as Lagrangian 
averaging or balanced models. We also briefly mentioned the interesting case of time-dependent 
background profiles, which was discussed in [DurOS]. These models are useful for driving local 
area models using profiles from global model output, or to model seasonal variations or regional 
climate change. In this case, potential vorticity will be conserved but we expect changes in domain- 
integrated energy in general, since there is now an explicit time- dependence in the Lagrangian. It 
is interesting to consider the case where the time-dependence in the background flow is slow, since 
in that case we expect the energy to be almost invariant; we will study this problem in future 
work. 

In the second half of the paper, we adapted this derivation to obtain vertical slice soundproof 
models. These vertical slice models are used to model front formation, and allow for full three 
dimensional velocity vectors, but with all fields being ^/-dependent except for potential tempera- 
ture which is assumed to have a constant gradient in the ?/-direction. Unlike the incompressible 
Boussinesq model, our slice models do not recover solutions of the full 3D equations since they 
have linearised the equation of state under the assumption of small 6 variation in the ^/-direction. 
However, we believe that they can play a role in the validation of numerical schemes and closures 
as part of the programme of asymptotic limit analysis described in [Ciil07], in which numerical 
solutions are compared with numerical solutions of the limiting semi-geostrophic model, which are 
most easily obtained for models in slice configuration using a geometric algorithm. Initial work on 
this programme for Boussinesq models is shown in [VCC13]. A requirement for such tests is that 
the model being tested has a conserved energy and potential vorticity. Both conservation laws are 
provided by these models, which can be obtained from minor adaptations of existing soundproof 
codes. 
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